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Abstract.  We  use  a  fully  self-consistent  three-dimensional  quantum  mechanical  transport  formalism  to  examine 
the  performance  of  InAs  based  quantum  wire  transistors  both  in  the  ballistic  limit  and  with  phonon  scattering 
included.  We  present  a  method  for  the  inclusion  of  polar  optical  phonon  scattering  as  a  real-space  self-energy 
tenn.  We  find  that  the  ballistic  perfonnance  of  the  devices  can  be  recovered  if  the  dopants  in  the  system  are  kept 
away  from  the  channel  entrance  and  exit.  When  dopants  are  present  at  these  key  points,  we  find  that  the  altered 
carrier  energy,  particularly  in  the  source,  has  a  significant  impact  on  the  device.  This  ballistic  recovery  is  aided 
by  the  fact  that  at  higher  energies,  polar  optical  phonon  scattering  loses  its  non-locality  which  leads  to  a  reduced 
scattering  rate  in  these  confined  systems. 


1.  Introduction 

It  has  been  a  well  established  fact  that  the  semiconductor  industry  is  aggressively  seeking  alternative 
device  solutions  in  order  to  continue  the  relentless  scaling  of  device  density.  One  of  the  most 
promising  of  these  new  technologies  is  the  tri-gate  quantum  wire  MOSFET  (Metal  Oxide 
Semiconductor  Field  Effect  Transistor).  Tri-gate  devices  offer  greatly  improved  electrostatic  control 
over  the  channel  of  the  device  when  compared  to  the  traditional  MOSFET  due  to  the  presence  of  the 
additional  gates  [1], 

Within  the  last  few  years,  there  has  been  a  renewed  interest  in  the  use  of  III-V  materials  instead  of 
silicon  for  future  generation  transistors  [2].  Previous  studies  have  shown  that  care  must  be  taken  in  the 
fabrication  of  these  devices  or  large  increases  in  the  output  drain  current  (kink  effect)  occur  [3-4], 
However,  it  has  been  shown  that  through  the  use  of  quantum  wires  with  discrete  subbands  for  the 
active  regions  of  the  device,  the  kink  effect  may  be  eliminated. 

In  this  paper,  we  present  the  results  of  a  three-dimensional,  self-consistent,  quantum  mechanical 
device  simulation  of  several  tri-gate  InAs  quantum  wire  MOSFETs.  Here  we  present  results  of  10  nm 
channel  length  devices  both  in  the  ballistic  limit  and  with  the  effects  of  impurity,  acoustic  deformation 
potential,  intervalley  (T— >X  and  F— >L)  and,  polar  optical  phonon  (POP)  scattering  as  separable  real 
space  self-energy  terms. 

2.  Device  Structure 

In  fig.  1,  we  display  a  schematic  of  the  device  geometry  for  an  InAs  MOSFET  with  a  9.69  nm  channel 
taken  to  lie  in  the  x-y  plane  (the  z-axis  is  normal  to  the  plane  shown).  The  exact  device  dimensions 
(multiples  of  the  lattice  constant)  have  been  included  in  this  simulation  to  aid  in  the  inclusion  of  the 
discrete  dopants.  The  thickness  of  the  InAs  layer  is  9.09  nm.  The  source  and  drain  of  the  device  are 
n-type  with  a  doping  density  of  6  x  101S  cm"3,  while  the  channel  of  the  device  is  considered  to  be  p- 
type,  but  left  undoped.  The  gate  material  is  assumed  to  be  platinum  and  the  gate  oxide  on  each  side  is 
composed  of  1  nm  of  hafnium  oxide  (HICh).  Underneath  the  device,  we  have  assumed  a  generic 
insulating  substrate.  Once  the  device  geometry  is  defined,  the  InAs  lattice  is  scanned  and  the  dopants 
are  distributed  according  to  the  method  presented  in  [5].  Following  the  distribution  of  the  dopants, 
they  are  then  mapped  back  onto  the  grid  of  the  simulation  mesh  and  the  initial  self-consistent  Poisson 
solution  is  obtained.  In  this  case,  the  full  Coulomb  potential  of  the  dopants  is  incorporated.  Then,  the 
solution  of  Poisson's  equation  for  the  local  potential  is  no  longer  smoothly  varying  in  the  source  and 
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the  drain  of  the  device.  The  inclusion  of  discrete  dopants  causes  the  fonnation  of  potential  variation  in 
the  source  and  drain.  The  density  throughout  the  device  is  calculated  using  a  variant  of  the  recursive 
scattering  matrix  method  which  solves  the  Schrodinger  equation  in  the  effective  mass  approximation 
as  described  in  [6].  In  short,  this  method  provides  us  with  a  three-dimensional,  fully  quantum 
mechanical  method.  Since  the  transport  calculation  is  performed  in  real  space,  the  different  excited 
modes  in  the  system  are  automatically  coupled  unlike  in  other  simulations  [7].  In  order  to  achieve  self- 
consistency,  the  density  obtained  from  the  transport  calculation  is  then  updated  using  Broyden’s 
method  [8]  and  a  new  guess  for  the  potential  is  obtained  through  the  solution  of  Poisson’s  equation. 
The  process  is  repeated  until  a  desired  level  of  convergence  is  obtained.  It  should  be  noted  that  all  of 
the  simulations  presented  here  have  been  performed  at  300  K. 
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Figure  1 :  Schematic  of  the  device  under  consideration  taken  in  the  xy  plane. 


3.  Polar  Optical  Phonon  Scattering 

The  inclusion  of  separable  scattering  mechanisms  is  not  a  new  concept.  In  previous  work,  acoustic 
and  optical  phonon  processes  have  been  derived  [9]  and  so  the  work  will  not  be  repeated  here. 
However,  we  are  now  interested  in  including  polar  optical  phonons  which,  to  this  point,  have  not  been 
treated.  We  begin  by  noting  that  to  use  the  method  presented  in  ref.  9,  we  are  assuming  that  the 
scattering  is  weak  relative  to  the  energies  in  the  system,  and  that  we  are  only  interested  in  the  steady 
state  limit.  Therefore,  we  can  use  the  Fermi  golden  rule  expression,  equivalent  to  a  first-order,  non¬ 
self-consistent  Born  approximation,  for  each  of  the  scattering  processes  and  generate  a  real  space  self¬ 
energy  from  it.  In  particular,  we  note  that  the  imaginary  part  of  the  self-energy  term  is  related  to  the 
scattering  rate  and  it  is  the  latter  scattering  rate  that  we  wish  to  calculate  [9].  This  will  result  in  an  x 
directed  momentum  which  is  related  to  the  carrier  energy  in  the  quantum  wire.  Since  many  of  the 
details  of  the  derivation  of  similar  scattering  rates  are  outlined  in  ref.  9,  we  will  show  only  the 
important  changes  in  the  derivation.  In  this  paper,  we  use  the  Fermi  golden  rule  for  scattering.  This  is 
treated  in  many  other  places  and  the  only  modification  that  we  need  is  to  account  for  the  transverse 
modes  of  the  quantum  wire.  Therefore,  we  begin  with  the  general  form 
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The  delta  functions  in  eq.  (1)  serve  to  conserve  the  energy  in  the  process  of  the  interactions  of  the 
carriers  with  the  polar  optical  phonons,  ip f  -(y,z)  is  the  transverse  two-dimensional  wave  function  in 
the  particular  slice  under  consideration.  We  assume  in  this  treatment  that  we  are  dealing  with 
parabolic  bandstructure.  Taking  a  closer  look  at  the  terms  contained  within  the  delta  functions,  we 
realize  that  based  on  this  assumption 
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In  eq.  (2),  and  Eff  are  the  initial  and  final  energies  at  the  bottom  of  their  respective  bands.  We 
begin  to  simplify  eq.  (1),  by  first  examining  the  longitudinal  terms 
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We  now  expand  eq.  (3)  using  the  results  from  eq.  (2) 
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This  sum  maybe  simplified  further  by  taking  the  Fourier  transform, 
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Representing  eq.  (5)  as  a  summation,  we  have 
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At  this  point  we  examine  the  relationship  for  E: 
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Eq.  (8)  is  quadratic  in  qx  and  may  then  be  solved  to  yield  a  solution  for  qx  wliich  is  then  substituted 
back  into  eq.  (7)  making  the  term  in  the  absolute  value  independent  of  qx  as  shown  below  regardless  if 
the  phonon  process  under  consideration  is  emission  or  absorption.  Therefore,  we  may  group  all  of  the 
variables  that  are  independent  of  qx  into  a  simple  dummy  variable  a.  With  this  pairing,  we  may  reduce 
our  summation  in  eq.  (7)  to  reach 
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Eq.  (9)  is  most  easily  evaluated  by  making  a  transformation  to  cylindrical  coordinates  to  obtain 
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We  evaluate  the  integral  in  eq.  (10)  using  contour  integrations  to  arrive  at  a  final,  simplified 
expression  for  the  sum  over  the  longitudinal  wavevectors 
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The  expression  in  eq.  (11)  is  substituted  back  into  eq.  (1)  to  find  an  overlap  integral  that  is  similar  to 
both  the  one  found  in  ref.  9  and  in  most  textbooks 
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Where  in  eq.  (12)  we  have  used  the  following  definition, 
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Following  the  usual  procedure,  we  replace  the  summation  over  the  final  momentum  states  with  an 
integration 
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00  00 

2^2 L-\dk'  =  \?w(E')dE' 

k'  —oo  0 


1 


^  polar  j . 


*  T  2 

mxLe 

3  27t3yco0  l'J 


(14) 


(15) 


However,  we  are  still  in  the  mode  space  representation  and  the  longitudinal  momentum  is  not  yet  a 
viable  operator.  To  remedy  this,  we  take  the  inverse  transform.  This  essentially  involves  solving  a 
contour  integral  for  both  the  emission  and  absorption  cases  which  results  from  eq.  (15).  This  gives  us 
our  final,  simplified  and  summarized  results  for  absorption  and  emission  of  POPs, 
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This  is  a  rather  expected  result  as  we  now  see  the  manifestation  of  the  non-local  nature  of  POP 
scattering  as  we  see  a  clear  dependence  on  the  grid  in  eq.  (16).  To  utilize  this  in  the  transport 
calculation,  we  must  use  the  following  unitary  transformation  to  convert  this  form  to  the  site 
representation  for  inclusion  in  the  Hamiltonian, 
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where  U  is  a  unitary  mode-to-site  transformation  matrix.  The  unitary  matrix  U~  results  from  the 
eigenvalue  solutions  in  the  transverse  slice  and  are  composed  of  the  various  eigenfunctions  in  the  site 
basis.  Hence,  it  represents  a  mode-to-slice  transformation. 


4.  Results 

In  fig.  2  we  plot  the  /,/- V,,  curves  for  the  ballistic  and  quasi-ballistic  case  for  two  different  dopant 
distributions.  In  fig.  2(a)  we  can  clearly  see  that  the  inclusion  of  the  phonon  processes  has  a  profound 
effect  on  the  device  performance.  This  degradation  in  performance  can  be  directly  linked  to  the 
inclusion  of  the  POP  scattering.  The  intervalley  transitions  (T— >X  and  T— >L)  are  quite  limited,  in  our 
case,  as  these  satellite  valleys  are  offset  by  0.73  eV  and  1.02  eV  respectively.  Therefore,  only  at  high 
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voltages  would  the  carriers  be  capable  of  reaching  these  valleys.  Interactions  with  the  acoustic 
phonons  are  still  present,  but  never  quite  reach  large  enough  energies  to  effect  the  transport  in  a 
significant  manner.  While  it  is  rather  expected  that  POP  scattering  causes  significant  reductions  in 
device  performance,  the  role  that  quantum  interference  appears  to  have  in  device  operation  is 
unexpected.  In  fig.  2(b),  we  plot  the  I,r  Vg  curves  for  the  ballistic  and  quasi-ballistic  case  for  a  different 
dopant  distribution.  Here  we  find  that  the  device  nearly  recovers  its  ballistic  device  performance. 
This  ballistic  recovery  is  attributed  to  two  factors,  localization  of  the  polar  optical  scattering  process 
and  differences  in  dopant  locations.  As  the  energies  present  in  the  system  increase,  we  find  that  the 
POP  scattering  becomes  relatively  localized  which  leads  to  reductions  in  both  the  absorptive  and 
emissive  coupling.  Therefore,  higher  energy  modes  are  subjected  to  reduced  scattering.  Perhaps  more 
interestingly,  we  find  that  quantum  interference,  of  the  carriers  with  the  dopants,  is  not  quenched  in  the 
presence  of  scattering.  The  device  in  fig.  2(a)  has  dopants  present  close  to  the  source-channel  and 
channel-drain  areas  which  significantly  impacts  carrier  energies  of  the  electrons  entering  the  channel. 
The  device  in  fig.  2(b)  does  not  have  dopants  near  the  channel  entrance  and  exit  areas,  and  therefore, 
has  reduced  quantum  interference  effects. 


(a)  (b) 


Gate  Voltage  (V)  Gate  Voltage  (V) 

Figure  2:  Id-Vg  curves  under  ballistic  (solid)  and  quasi-ballistic  (dotted)  conditions  for  an  identical 
dopant  distribution  where  (a)  dopants  are  present  close  to  the  channel  interfaces,  and  (b)  dopants  are 
removed  from  the  channel  interface,  respectively. 

5.  Results 

In  this  paper,  we  have  presented  results  of  self-consistent  three-dimensional  fully  quantum  mechanical 
simulations  of  10  nm  InAs  tri-gate  quantum  wire  transistors.  We  have  additionally  presented  a  method 
for  the  inclusion  of  POP  as  a  real-space  self-energy  tenn.  We  find  that  the  ballistic  performance  of  the 
devices  can  be  recovered  if  the  dopants  in  the  system  are  kept  away  from  the  channel  entrance  and 
exit.  When  dopants  are  present  at  these  key  points,  we  find  that  the  altered  carrier  energy,  particularly 
in  the  source,  has  a  significant  impact  on  the  device.  This  ballistic  recovery  is  aided  by  the  fact  that  at 
higher  energies,  POP  scattering  loses  its  non-locality  which  leads  to  a  reduced  scattering  rate  in  these 
confined  systems. 
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